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Abstract. Starting from microscopic interaction rules we derive kinetic mod- 
els of Fokker-Planck type for vehicular traffic flow. The derivation is based on 
taking a suitable asymptotic limit of the corresponding Boltzmann model. As 
particular cases, the derived models comprise existing models. New Fokker- 
Planck models are also given and their differences to existing models are high- 
lighted. Finally, we report on numerical experiments. 

1. Introduction 

Kinetic models for traffic flow provide a description of vehicular traffic based 
on spatial, temporal and velocity information and have been subject to intense 
discussion in recent literature, see e.g., [H El [9j \TU\ and [HJ H2]- In this paper 
we contribute to this discussion in the following way: we introduce microscopic 
interaction rules and derive existing and new equations using formal asymptotics 
and limits. This approach provides insight in the dominant terms of kinetic models 
and shows possible extensions. 

Microscopic interaction rules yielding Boltzmann type traffic flow equations have 
been analyzed in [2]. Therein, the models by Klar et. al. [TJ El El HO] could be derived 
from basic rules. However, the important class of Fokker-Planck type traffic flow 
models as introduced by Illner et. al. [11 could not be obtained. The present 
work now closes the remaining gap. We introduce general microscopic interaction 
rules as well as a derivation leading to the Fokker-Planck type models of Illner et. 
al. Further models corresponding to modified interaction rules are also discussed. 
In this work we do not want to discuss further the validity of such Fokker-Planck 
models (see [TT] for a detailed discussion). We only remark that these models 
are particularly interesting from the computational viewpoint since they keep the 
kinetic information on vehicles without requiring the evaluation of the expensive 
velocity integrals as in Boltzmann-like models. 

The main target is the single-lane, spatial homogenous Fokker-Planck type equa- 
tions for traffic flow. The prototype example introduced ad-hoc in [TT] is given by 

(1) d t f + d v (B[f]f - D[f]d v f) = 0, 
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where f(v,t) is the density of vehicles with velocity v at time t. Its motivation is as 
follows: B is a hcuristically defined braking and acceleration operator depending on 
moments of /, i.e., macroscopic quantities, like number density p or mean velocity 
u. They are defined as 



The underlying assumption on both operators is that drivers only observe averaged 
quantities. The diffusive term is introduced to model the fact that drivers may 
not observe the averaged moments correctly. Furthermore, this operators allows for 
trivial equilibria corresponding to synchronized flow. Examples for the operators 
B[f] and D[f] are given in [TT] and [3] 



Herein , c\ are positive constants corresponding to braking and acceleration, \ 
is the indicator function. For convenience we normalize density and velocity such 
that pmax — u max — 1. The function a(p,u) ensures that the boundary conditions 
in the stationary case are met. It vanishes for large and small values of p and u. An 
explicit form can be found e.g. in [llj . It can be changed to a(p, u)\v — m| 7 + £, where 
e is a residual diffusion to prevent discontinuous equilibria, see [3]. Furthermore, 
the function q(p, u) is motivated by additional lane-changing activity. It is modeled 
as 



The coefficients 7 and 6 satisfy 0<7 + <5<3in order to obtain well-defined steady 
states. 

In this paper, starting from microscopic considerations, we obtain a microscopic 
motivation for the given form of braking/ acceleration and diffusion operators in 
the Fokker-Planck model. The microscopic rules will also lead to different but still 
qualitatively similar Fokker-Planck type models. Hence, in the sequel we start from 
a microscopic model with precise interaction rules between corresponding cars and 
derive on this grounds kinetic models of Boltzmann type. A suitable asymptotic 
limit then permits to recover the desired Fokker-Planck type models. We give the 
presentation for a general interaction of rule and state the corresponding equations 
for several well-known examples. Finally we report on numerical results for the 
different Fokker-Planck models. 

2. Microscopic interaction rules and Boltzmann models 

The interaction rules in traffic flow are different from the usual gas-dynamics 
case. In particular, traffic "particles" are anisotropic, since a car that drives be- 
hind another one and speed up does not necessarily make the person in front go 
faster or slower. Additionally, acceleration and braking scenarios are not necessarily 
symmetric to each other. 

We will present the derivation of the kinetic models for a general microscopic 
interaction rule. The leading car is assumed to have a pre-interaction velocity v 
and a post-interaction velocity v'. We distinguish two different scenarios: if the 
leading car is slower than the one in front of it the leading car accelerates, cf. 



P 




(2) 




d v {(a(p,u)\v - u[')d v f}. 
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(l3al) . On the other hand, if it is faster it will decelerate or brake, cf. ([3b)) . A 
simple model is as follows: we compare our velocity v with the velocity w of the 
car in front. If we are faster, we brake according to the relative velocity, if we are 
slower we accelerate. Since we cannot precisely estimate the leading car's velocity 
and we will made mistakes when braking and accelerating. The modeling includes 
therefore some additional parameters and a noise term proportional to braking and 
acceleration intensity characterized by a random variable £ . 
The general form of the microscopic interaction rule is 

(3a) Acceleration: v = v + A(V A - v) + £v(v)(Va ~ v) K , v < W 

(3b) Braking: v' = v - B(v - V B ) + fy(v)(v - V B ) R , v>W 

where W G (0, 1) represents the estimate on the leading car velocity w made by the 
driver. The coefficients Va £ (v, 1], Vb £ [0, v) represent the desired velocities when 
drivers accelerate and brake. Thus the natural choice is W = W(w) with W(-) a 
monotone non decreasing function of in. Beside the obvious choice W = w other 
choices for W, although less realistic, are interesting since they lead to simplified set 
of equations. For example W = u where drivers are able to estimate only the mean 
velocity of the flux. Similarly Va = Va{v, W) and Vb — Vb{v, W) may depend on 
the actual velocity v and the estimated relative velocity W. 

To complete the interaction we leave the velocity of w of the car in front un- 
changed, i.e., 

w = w. 

In ([3]) A and B are non negative constants in (0, 1), k > 1 is a power that 
calibrates the dependence of the noise from the braking/acceleration dynamic, < 
v{v) < 1 is a function that takes into account that noise is also proportional to 
the actual driver speed and £ is a random variable with zero mean and variance a 2 
distributed accordingly to 

The above microscopic relations include as particular cases several well-known 
models. We give some examples. 

• Example 1: The model by Illner, Klar & al. [6]. In this model we assume 
that individual estimates correctly the leading car's velocity and that ac- 
celeration and braking are bounded by the maximal and minimal velocity 
instead of the velocity of the leading car 

W = w, V A = 1, V B = 0, v(v) = 1, k = 1. 

• Example 2: A generalized model by Giinther, Klar & al.[4]. Drivers in- 
teract via symmetric braking and acceleration scenarios governed by the 
leading car and can estimate correctly the leading car velocity 

W = w, Va = w, Vb = w. 

Note that in the original model no diffusion was present v{v) = 0. 

• Example 3: If we assume that individuals react with respect to the aver- 
age velocity and that acceleration and braking are bounded by the average 
velocity we have the set of parameters 

W = u, Va — u, Vb = u. 
This model as we will see permits to derive the Fokker-Planck equation ((2|) 
from the microscopic dynamics in a suitable asymptotic scaling. 
Some remarks are in order. 



4 



MICHAEL HERTY AND LORENZO PARESCHI 



Remark 1. 

(1) In the microscopic model it is possible to add a further term independent 
of the other drivers behavior 

(4) v' = v + E(4>{v)-v). 

From a modeling point of view this corresponds to free flow traffic ad- 
justing towards a certain desired velocity distribution tp(-). To keep the 
derivation simple we will ignore the presence of this free flow term in the 
sequel. 

(2) Note that in ([3]) we require that the post-interaction velocities satisfy v' S 
(0,1). We can show that this is guaranteed by our assumptions on the 
parameters. Let us show this for the acceleration term. Since < Va — v < 1 
we have 

v' < v + A{V A -v) + \i\v{v){V A -v) 
= v (l-A-\Z\»(v)) + (A + \Z\v(v))V A 
< V A 

if |£| < 1 — A. Note that for A close to 1 the domain of £ vanishes. One can 
show that independently of A we have v' < 1 if |£| < (1 — Va)/Va- However 
in the sequel we will be interested in the behavior of the model for small 
values of A and B and for such reason we will take simply |£| < 1 — A. 

It is easy to see that v' > v if |£| < A. Note however that in this case for 
small values of A the domain of £ vanishes. If we assume that v(v) < v we 
get 

v' > v + A(Va-v)-\^\u{v){V a -v) 

> v + A{V A -v)-\£\v 
= v{l- A - |£|) + AV A 

> AV A 

if |£| < 1 — A which permits to pass to the limit A — > keeping £ finite. 

Similarly for the braking term we have v' € (Vb, 1) provided that |£| < 
1 — B and v(v) < 1 — v. Thus if we want to preserve the lateral bounds 
of v. To this end the presence of the function v(v) is of importance, v is 
such that it vanishes at the extreme values of the velocity v = and v = 1. 
This can be achieved by any function such that v(v) < min{u, 1 — v} for 
v G (0, 1). For example we might take 

!/(«) = i(i-|2«-l|) 

or a quadratic function like 

v{v) = v(l — v). 

This implies that the diffusion term in (I3a[) - (l3b[) vanishes in correspondence 
of the desired speeds Va,b and at the extremal values and 1 of the velocity 
domain. 

As mentioned in the introduction a general Boltzmann-type formulation for ve- 
hicular traffic was derived in 0. For simplicity we restrict to the case W = w in 
([3]) we will come back to other choices at the end of the next section. 
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In weak form we can write for the car distribution f(v,t) the kinetic equation 

/•OO 

(5) d t / }{v)4>{v)dv = Q A + Q B 

Jo 

where 

(6a) Qa := f 1 / V^^^/W/HW''')-*))^^^ 
Jo Jv Jr 

(6b) Qb := f [ V [ ^^^f^Kwmv^-^d^dwdv. 

Jo JO JR 

Additionally, as in [2] we assume that the function (3 is given by 

g A,B = f v(0\ v - w\QaX[o,i](v') v < w 

P(v,w)^{v>,w>) \ v(0\v - w\qBX[o,i]{v') V > w, 

where \ is the indicator function. In the general case © the collision kernel has a 
weight of \v — W\. 

From (f6a[) - (|6b[) conservation of the total number of vehicles is obtained for <fi(v) = 
1 which represents the only conservation property satisfied by the equation. If we 
consider the case 4>{v) — v (this corresponds to the behavior of the mean velocity) 
we obtain 

|>00 f 1 f 1 {• 

dt / f(v)v dv = A / / Pf v , w) ^ iv/ , w/) mfH(VA - v) dw dv + 



Jv 



P( V , W )^(v>, w i)mfM&(v)(v A -vydtdwdv 

U I r J/(v^(v'^)f(v)f(w)(v - V B ) d^dwdv - 
B tr„.\ tr...\t:..f..\r«. ir ^fc 



^0 



P(v^(v. w ')f(v)f(w)t;v(v)(v - V B ) K dtdwdv. 

10 JO JR 

Hereby we assume that the random variables in the acceleration and the braking 
terms take values on the sets (A — 1, 1 — A) and (B — 1, 1 — B) respectively. Now 
the second and the fourth integrals vanish since £ has zero mean and so we have 

roc rl /"l 

®t / f(v)vdv = A / / \v - w\f(v)f(w)(VA — v) dw dv — 

JO JO Jv 

pi pv 

B I j \v — w\f(v)f(w)(v — Vb) dw dv, 
Jo Jo 

which represents the change in mean velocity due to acceleration and braking. 
If we assume A = B = 7 and Va = Vg = V we can write 

/>oo />1 pi 

®t / f(v)vdv = j / / \v — w\f(v)f(w)(V — v) dw dv. 
Jo Jo Jo 

In Example 2 we have Va = V B = w. In this case a straightforward computation 
shows that the mean velocity u is conserved. This situation has much in common 
with a granular gas dynamic where 7 is related to the coefficient of restitution, and 
the large time behavior of the system in absence of diffusion {v(v) — 0) is described 
by synchronized traffic given by a Dirac delta f(v)=5(v — u) where all cars have 
the same speed u. Note that the synchronized traffic state is an equilibrium state 
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of the system even in the general case when A ^ B and v (v) ^ since the noise 
term vanishes when all cars have the same speed. 

One of the major drawbacks of Boltzmann-type models is their excessive com- 
plexity which makes it difficult the development of numerical simulation methods 
as well of analytical tools for the computation of steady states. 

In order to obtain simpler kinetic Fokker-Planck type equations we use a suitable 
asymptotic limit similar to the technique presented in |13] . In the latter presentation 
only a global and symmetric interaction rules for a Maxwellian-like kernel have been 
considered. Other similar interaction rules in the context of opinion formation 
has been studied in [2]. We emphasize that in the following, the interaction is 
anisotropic and the collision kernel is of hard-spheres type. Therefore it is different 
to [13l [14] and also to standard approaches in gas-dynamics. 



3. Fokker-Planck type asymptotic models 

First let us recall some definitions. Wc define Aio the space of all probability 
measures in R + and by 

(7) M p = 1 77 6 Mo : jf \tf V (0 d£ < +oo,p > oj , 

the space of all Borel probability measures of finite momentum of order p, equipped 
with the topology of the weak convergence of the measures. Let .7p(R + ), p > 1 
be the class of all real functions on R + such that g(0) = g'(0) = 0, and g^ m \v) is 
Holder continuous of order 5, 

g( m > s = sup — \ y < oo, 

the integer m and the number < 5 < 1 are such that m + 5 = p, and </ m ) denotes 
the m-th derivative of g. 

In the rest of the paper we will assume that the symmetric probability density 
fy(£) which characterizes the transition rate belongs to A^2+a, for some a > 0. 
Moreover, to simplify computations, we assume that this density is obtained from a 
given random variable Y with zero mean and unit variance, that belongs to M.2+a- 
Thus, rj of variance a 2 is the density of aY. By this assumption, we can easily 
obtain the dependence on a of the moments of 77. In fact, for any p > 2 such that 
the p-th moment of Y exists, 

\t\*T,(Qdt = E(\*Y\*)=o*E(\Y\>). 

To skip inessential difficulties, that do not change the following analysis, we 
suppose that the random variables in the acceleration and the braking terms take 
values on the sets (A — 1, 1 — A) and [B — 1, 1 — B) respectively. We also assume 
that (f> £ J" 2+5 (K + ) for a > 6 > 0. 

In order to derive a Fokker-Planck equation from ([6]) we need to consider an 
appropriate scaling of the equations. We set for e > 

(8) r = et, f(v,r)=f(v,t). 
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and assume that there exists non-negative constants ca,cb and cd with the fol- 
lowing properties 

A B „ a 2 

(9) hm — = ca, hm — = cb, hm — = co- 

A^0 : £^0 £ B->0,e^0 £ " cr->0, e->0 e 

These assumptions corresponds to an asymptotic situation where each interaction 
produces very small velocity variations and in the limit we preserve the main features 
(deterministic acceleration/braking and random noise) of the dynamic. 
The scaled kinetic equation reads 

(10) d T f f(v)<j>(v)dv = -(Q A (f,4,,r ] )+Q B (f^,r l )). 

Jo e 

For small values of A, B and a 2 we have v' « v and so we can consider a Taylor 
expansion in the collision kernels Qa and Qb of <j)(v') around v. Let us consider 
here the acceleration term (l3al) 



cp(v') - 0(«) = (A(V A -v) + tu(v)(V A ~ v) K )Mv) + 
\{A{V A -v)+ fr(v)(V A - v) K ) 2 <t> vv (v), 
where for some 9 S (0, 1) we have 

i = (1 -6)v + 0v'. 

Now we insert this expansion into the acceleration operator Qa and compute the 
limit for e — > 0. Since rj(£) has mean zero and vanishing variance we get 

lim -QaUAiV) = QaIca [ [ (w — v)(Va —v)f{v)f(w)<p v (v)dwdv. 
e L Jo Jv 

c D / (w~v)v(v) 2 {V A -v) 2K f(v)f(w)(/) vv {v)dwdv 



Jv 



where 



2 



lim -R(f, farj) 
e->o e 



l r i 



R{fA,rj) = \q A c D I I I r](0(w-v)v(v)(A(VA-v)+&(v)(VA-vrY- 
JO Jv Jr 

■f{v)f{w){(f)vv{v) - 4>w{v)) d^dwdv. 
Since </> £ J r 2+s('^+), and \v — v\ = 6\v' — v\ 

(11) \<t> vv (v) ~ <M«)I < \\<t>vvh\v - vf < Uvvhlv' - V\ S . 

Hence 

\R(f,c/>,r,)\ < qAC D 1 ^^ f f ( - vHv) ■ 

z Jo Jv JR 

■\A(V A -v) +£u(v)(Va - vf\ 2+s f{w)f{v)d^dwdv. 

By virtue of the inequality 

\A(V A -v) + iv{v){V A - vT\ 2+s < C 5 (A 2+S + \Z\ 2+S ), 

and using the fact that v is compactly supported in (0, 1) we finally obtain the 
bound 

(12) \R(f,cb, v )\ < qACD M^k Cs p 2 (A 2 + s + [r,m\ 2+5 dC 
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Since 77 is a probability density with zero mean and a 2 variance, and rj belongs to 
M2+a, for a > 5, 

[ l£| 2+ VCR = E (\aYf +S ) = o^E (|y| 2+5 ) , 

and E (\Y\^ +S ^j is bounded. Using this equality one shows that jR(f,<f),r)) con- 
verges to zero as e — > 0. 
Thus we obtain 

lim -Q A (f,(f),ri) = qA^CA ( / (w - v)(V A - v)f(v)f(w)<f> v (v) dw dv. 

6 I JO Jv 



1 r 1 * 



2 



Jv 



+ {w-v)v{vY{V A -vY K f{v)f{w)^ vv {v)dwdv 



By similar arguments for Qb we get 

lim -Q B (/,</>, rj) = q B \ c B [ { (v - w)(V B - v)f{v)f{w)<j) v {v) dw dv. 
e I Jo Jo 

+ I c dJ o I {v~w)v{vf{V B -vf K f(v)f{w)^ vv (v)dwdvy 

In order to present the strong formulation of the limiting kinetic equation we intro- 
duce the following constants 

(13) ca = c A q A , Xa = c D q A , X B = c D q B , c B = c B qB- 

Summarizing, the above computations we observe that in the limit e — > we obtain 
the following Fokkcr-Planck equation given in strong form by 

d T f + d v < f I c A / f(w)(w - v)(Va - v)dw - c B f(w)(v - w)(v - V B )di 
(14} V Jv J ° 

^d vv ^u(vff (x A J f(w)(w - v)(V A - vf K dw + X B J f(w)(v - w)(v - V B f K dv< 
which can be rewritten in more compact notation as 

(15) d T f + d v (f(v) f \v-w\L(v,w)f(w)dw 

(16) V Jo 

^d vv (v(v) 2 f(v) J \v - w\D(v,w)f(w)dw ) 



where 



c A (V A -v) v<w, ( X a \V a -v\ 2k v<w 



1 c B {V B -v) v>w, v ' ; 1 X b \V b -v\ 2k v>w. 
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Furthermore, it is useful to introduce the partial moments for fast and slow cars 
as follows 



Ps = f(w)dw, pf = f{w)dw, 

JO Jv 



1 /"" . . 1 



1 



ug — — / f(w)wdw, up = — / f(w)wdw, 

PS Jo ' PF J v 

i r if 1 

T s = — f(w)(v -w)(v -V B )dw, T F = — f(w)(w-v)(V A -v)dw, 

PS Jo PF J v 

w£ = —[ f(w)(v-w)(v-V B ) 2K dw, W$ = — [ f{w)(w-v){V A -vf K dw. 

PS Jo PF Jv 

In this way we can rewrite the limiting Fokker-Planck equation in the form 
(17) d T f + d v (/ (c A p F T F - CBPsTs)) = ^d vv (v{v) 2 f (\ APF W£ + \ B p S W§)) 

For the interaction rules presented in the previous section we have the following 
equations in the limit e — > 0. 

• Example 1: For the model by Illner, Klar & al. we get 



w)dw 



(18) 
1 
2 



d T f + d v j/ ^ca(1 - v) J f(w)(w -v)dw - c B v J f(w){v- 

A^(l — v ) 2 j f{w)(w — v)dw + \BV 2 J f(w)(v — w)dvj S j | 



\d vv { v(v) 2 f 



or using the partial moments 

d T f + d v (/ (ca(1 - v)(p F u F - p F v) - c B v{p s v - psus))) = 
^d vv {v{v) 2 f (Aa(1 - v) 2 (p F u F - p F v) + X B v 2 (p s v - psus))) 



• Example 2: For the generalized model by Giinther, Klar & al. we obtain 

(19) 



d T f + d v |/ (c A J f(w)(w - vfdw - ~c B £ f(w)(v - w) 2 dw^j J 
\dw ^{vf! [\a J f(w)(w - v) 2+k dw + X B jf 



f(w)(v -w) 2+k dw 

• Example 3: For the model where individuals react with respect to the 
average velocity we have 

(20a) d T f + d v (/ (c A (v - u) 2 px{ v<u } - c B (v - u) 2 px{ v > u })) = 

(20b) l -d vv ( V (v) 2 f {X A (u - v) 1+2K p X{v<u} +X b (v- u) 1+2K px { v >u} )) 
Remark 2. 

• The evaluation of the partial moments which appears in examples 1 and 2 
can be done efficiently at a numerical level at the same computational cost 
O(N) of the computation of mass and momentum. This is of paramount 
importance in applications of the the Fokker-Planck models to realistic sim- 
ulations. 
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• The general Fokker-Planck equation (fT4l) derived above preserve the es- 
sential properties of the microscopic relations it came from. For example 
mean velocity is preserved only when Va = Vb = w and A = B. Simi- 
larly synchronized traffic is a possible equilibrium state of the system when 
Va = Vb = w. Moreover upper and lower bounds for the cars velocities are 
satisfied thanks to the presence of v(v). 

4. Relations with Illner, Klar, Materne model 

We analyze some of the derived models and compare them to Illner, Klar, 
Materne model. 

Let us at first consider Example 3 and assume that v(y) is independent of v and 

set 

X A = \b =■ A/^ 2 , 

i.e., q A = Qb- 

We can differentiate in the diffusion term and obtain 

drf + 8 V (f((c A (v - Uf + " v) 2 n P X{v<u } 

-(c B (v - uf + {l+ ^ K)X { v - u) 2k ) PX{v>u} )) = ~d v (d v f\u - v\ 1+2K \p) . 

By comparing this equation with the original model by Illner et. al. © we observe 
the following. Since we did not include lane changing in the microscopic interaction 
rules we can only compare with (J5]) in the case P = 0. The previous derivations 
hold true also in the case when the quantities ca, cb and A are p and u dependent. 
Hence, if they are chosen such that 

l + 2n x jl-p _ 1 + 2k j x . . , „ 

{®Ah o — A = c A , cb H — A = c B , \ = a(p,u)p, 1 + 2k = 7 

2 p 2 

we exactly recover the model by Illner et. al. © for 7 = 3. In the cases 7^3 
we can recover the exponent in the diffusion but obtain additional drift terms to a 
power of 7 — 1. We also observe that we immediately obtain the lower bound on 
7 > I for positive k. 

Moreover the diffusion term is degenerate when v — u. Adding a modified 
diffusion such as <j{p, u) + e will off course remove the degeneracy but destroys the 
possibility of synchronized traffic as discussed in [3] . 

Let us now consider Example 2. We approximate the partial moments in (|I9[) 
close to synchronized traffic 

f = S(v- u). 

At synchronized traffic the operators T5 and Ws collapse to 

Ts = ~(u - v) 2 X{ v >u}, T F = -(u- v) 2 X{ v <u}, 
p p 

W s = -(u - v) 1+2K X{v>u}, W f = -{u - v) 1+2K x{v<u}, 
p p 

and under the additional assumption A^i = =: A and using the same choices 
for the constants ca,cb and A as before, we again obtain ([2]) in the case P = 0. 
Therefore, the model (flU)) can be seen as a generalized equation comprising © and 
at the same time linking the coefficients appearing in diffusion and braking terms 
to microscopic interaction rules. 
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Note however that in the computation reported above to connect our Fokkcr- 
Planck models with Illner & al. model we assumed that v{v) is constant. Thus 
there is no mechanism to prevent the solution to violate the lateral bounds for the 
car velocity. 

We discuss some particular cases. The limit case in Illner's model for the diffusion 
coefficient is 7 = 1 which implies k — in the microscopic rules ([3]). This implies 
a constant diffusion in the microscopic interaction rules ([3]). This seems at least 
questionable. On the other hand taking k = I corresponds to 7 — 3 being the upper 
limit case for Illner' s model. In this case the noise is directly proportional to the 
difference in the velocities. 

Finally, we compare the diffusion terms of ([2]), i.e., 

d v (d v fa(p,u)\u - w| 7 ) 
and (TH)]) in the case A^ = A B = A, i.e., 

d v {dvf^Y^- f{w)\w - v\ 1+2K dv?j . 

We observe, that the Fokker-Planck derived from microscopic rules has a convo- 
lution integral in the diffusion term compared to the degenerated diffusion in the 
Illner model. Hence, in the steady states we expect a smoother distribution profile 
f(v). Furthermore, the degeneracy at v = u is not present in the second case. This 
is a remarkable feature of (fT9l) with respect to ([2]). 



5. Steady states and numerical experiments 

In this section we report on numerical results for the models ([15]). (|TT))) and 
(l20l) . In all models we set the cut-off in the diffusive part as v{v) — v(l — v) 
and hereby normalize the computational domain for v to (0, 1). We use N v = 100 
discretization points on (0, 1) for all computations. We apply the methods of lines 
with a time horizon of T = 250. The large time horizon is chosen in order to 
capture the steady state solution. When discretizing in v we use a simple first- 
order upwind discretization for the first-order derivatives and centered differences 
for the diffusion term. The final system of ordinary differential equations is then 
solved by an implicit Euler method. During simulations the total mass p — J fdv 
is preserved. The initial data for the simulation is chosen such that it vanishes at 
the boundary and given by /o = exp(— 25 (v — i) ). The initial mass is po « 0.35. 
For simplicity we set A = Xa = As in all simulations. 

At first we study the dependence of the time evolution of /o on the diffusion 
coefficient A if the dynamics is controlled by example 1 (If8l) . i.e., 



d t f + d v ( fc A (l - v) I f(w)(w - v)dw - c B v I f{w)(v - w)dw 



n 



^d vv ^{v) 2 fX^{l-v) 2 f{w)(w-v)dw + v 2 J f{w)(v-w)dw 

The braking and acceleration constants are cb = ca = I- We present plots of 
the isolines of the solution as well as a three-dimensional for diffusion coefficients 
between A = I and A = 50 in figure [T] We observe a concentration forming from 
the initial Gaussian curve. As we increase the diffusion parameter the center of the 
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concentration moves from v = s to the boundaries v = {0, 1}. A similar behavior 
is observed when simulating example 2, (|19p . i.e., 

d t f + d v (fcA J f(w)(w - v) 2 dw - c B J f(w)(v - w) 2 dv?j = 

\d m , (v(v) 2 f\J f(w)\w-v\ 2+K dwY 

Using the same parameters and setting k — we observe in figure [2] the same 
qualitative behavior. For small diffusion coefficients we observe a concentration at 
the center v = h and for larger coefficients at the boundary v = and v = 1. Next, 
we study the dependence of the time evolution on the other parameters present in 
example 2. In figure |3] we observe that the exponent k in (|19j) does not influence 
the qualitative behavior. We give a result for a moderate diffusion parameter and 
the case k = 1 and k = 2. However, the braking and acceleration constants have 
a strong influence on the position of the concentration. In the case k = 1 and 
A = 5 and |J > 1 we observe a shift of the point of concentration away from 
v = 4. Clearly, the obvious interpretation is: if the braking force is stronger than 
the acceleration force the equilibrium distribution must concentrate in a regime of 
lower velocity. The results are given in figure 2] and has to be compared with the 
corresponding part in the upper left of figure [T] where we used — = I . 




Figure 1. Example 1, numerical results for a diffusion coefficient 
of A = 1 (top left), A — 5 (top right), A = 10 (bottom left) and 
A = 50 (bottom right). 
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n 0.354383 0.354383 
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in 0.354383 0.354383 





J — 11 I 



Figure 2. Example 2, numerical results for a diffusion coefficient 
of A = f (top left), A = 5 (top right), A = 10 (middle left) and 
A = 25 (middle right) and A = 50 (bottom). 



Concerning Example 3, (|2"0"1) . a comparison between the derived model and the 
original model of Illner et. al. is given. Here, we compare the long term behavior in 
terms of the arising steady state. In case of Illner's model @ these can be computed 
explicitly [3] and are given by 

v > u 
V < u. 



(22) 



c exp 



c 2 exp 



-CbP 



(3 - j)a(p,u) 

-c^(l-p) 
(3 - i)a(p,u) 



(v — u) 
(u — v) 



3-7 



3-7 
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Figure 3. Example 2, numerical results for different interaction 
coefficients K = 1 (left) and k = 2 (right). 




0.1 0.2 0.3 0.4 0.5 0.G 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 4. Example 2, numerical results for cb/ca = 5 (left) and 
cb/ca = 10 (right). 



Here, it is assumed that P = 0. The constants c 1,2 are such that mass and momen- 
tum equation are satisfied, i.e., p = J fdv,pu = J fvdv. In case of Example 3 we 
also compute the steady states as 

„,._. 4i-((-^ ( "-"»^))-(-r^'"-»> 

c ( ( ^AP f_ n2+2k\ 1 I _ f °A 



ex P --^-(u-,)^ eX p /) -j^r^-t*)^ «<u 



^ 2 (w) r \\ 2 + 2k ' J J \ J v v 2 ^) 

Again, the constants c 1 ' 2 are such that p — f Q f(w)dw, pu — f Q wf(w)dw. The 
additional factor in the steady state is due to the fact that the diffusion is now 
velocity dependent. If the following relations hold true, both steady states coincide: 

(23) a(p, u) — pX, X — Xa — Xb, v(v) = 1, 1 + 2k = 7. 

6. Conclusions 

Starting from microscopic vehicles interactions we derived Fokker-Planck models 
for traffic flows as a suitable asymptotic limit of a Boltzmann dynamic. Among oth- 
ers the new Fokker-Planck models include the model recently introduced by Illner 
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& al. [2]. A remarkable feature of these models is to preserve synchronized traffic 
states and to avoid the degeneracy in the diffusion term when the vehicle velocity 
coincides with the mean velocity. These kinetic models are very promising for nu- 
merical simulation purposes since their computational cost is strongly reduced when 
compared to full Boltzmann models. Numerical results for the non homogeneous 
multilane case are under progress and will be presented elsewhere. 
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